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ABSTRACT 

High-energy 2-|im lasers have been incorporated in a breadboard coherent Doppler lidar to test component technologies 
and explore applications for remote sensing of the atmosphere. Design of the lidar is presented including aspects in the 
laser transmitter, receiver, photodetector, and signal processing. Sample data is presented on wind profiling and C0 2 
concentration measurements. 
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1. INTRODUCTION 


Recent advances in high-energy 2-pm solid-state laser technology have enabled lidar applications for long range wind 
sensing and C0 2 measurement. Wind profiling with 2 -pm lidar has already shown merit using low energy laser 
transmitters, and the addition of stronger energy lasers promises to boost the applicability of the coherent Doppler 
technique to farther distances and higher altitudes. 1 ' 2 ' 3 Research at NASA is aimed at developing lidar technology for 
satellite -based observation of wind on a global scale. The ability to profile tropospheric wind is a key measurement for 
understanding and predicting atmospheric dynamics and is a critical measurement for improving weather forecasting 
and climate model validation. 4 Homeland security applications have been proposed to provide wind data for accurate 
transport prediction of a plume release of a nuclear, biological, or chemical agent. 5 Safe and effective military 
operations are also dependent on weather forecasting and wind profiling. 

Another application for 2 -pm lidar is remote sensing of C0 2 by the differential absorption lidar (DIAL) technique. 
Atmospheric C0 2 concentration is of interest to climate change and global warming studies. The recent rate of increase 
of C0 2 is about 100 times higher than at any other time during the last 100 years. While the long-term trend of C0 2 in 
the atmosphere and its seasonal variations are well documented, the sources and sinks of C0 2 on continental and 
regional scales are not well understood. 6 A sink for C0 2 comprising more than 25% of the 7.1 petagrams (10 17 grams) 
of C0 2 that is released annually into the atmosphere, remains unaccounted for. 7 It is also known that both land and 
ocean sinks fluctuate with space and time. Knowledge of the spatial and temporal distributions of C0 2 and an 
understanding of the causes of these variations are essential in predicting future levels of C0 2 and their impact of 
climate change. Differential absorption lidar (DIAL) potentially offers several advantages over passive techniques for 
remote C0 2 sensing. These advantages include relatively high spatial resolution, high measurement specificity 
(avoiding interference from other gases), lack of dependence on external light sources, and relatively simple inversion 
methods compared to retrieval methods used in passive remote sensing. A challenge to lidar in this application is that 
the concentration measurement must be made with a very high level of precision and accuracy. To assist the 
understanding of many carbon cycle processes a C0 2 measurement accuracy of 1-2 ppm over 1-2 km vertical range 
intervals in the troposphere and column measurements with 1-2 ppm precision are ultimately needed. 



We report here on technology developments and lidar system demonstrations aimed at both wind and CO 2 applications. 
High-energy Ho:Tm:YLF and Ho:Tm:LuLiF laser transmitters are described along with aspects of injection seeding, 
wavelength tuning, and wavelength stabilization. Coherent heterodyne receivers are also discussed, including advances 
with an integrated photodiode/preamp in a dual-balanced configuration. A complete lidar system testbed has been built 
from these component and integrated with a real-time signal processing and display computer. Sample wind and CO 2 
measurements will be shown. Validar gets its name from the concept of “validation lidar,” in that it can serve as a 
calibration and validation source for future airborne and spaceborne lidar missions. Validar is housed within a mobile 
trailer for field measurements. 


2. LIDAR DESIGN 

The optical layout of Validar is shown in Figure 1 and a list of specifications is given in Table 1. Transmitter and 
receiver subsystems are each described separately, with the transmitter considered to consist of the pulsed laser, 
continuous wave (CW) laser, and the optics associated with coupling the two. The receiver is considered to be the 
transmit/receive switch, telescope, and photodetector. Not shown in Figure 1 is a scanner mounted through the roof of 
the trailer for steering the lidar beam. 



Figure 1: Optical layout of the lidar. 


2.1. Laser Transmitter 

The 2-Lim wavelength is advantageous in a coherent Doppler lidar for several reasons. First, it offers a high level of 
eyesafety with a maximum permissible exposure of 100 mJ/cm 2 . 8 Second, the 2-pm wavelength is strongly 


backscattered by aerosols in the atmosphere, from which the Doppler shift indicates wind speed. Third, holmium solid- 
state lasers can produce a high level of pulse energy at 2.05 pm. Last, C0 2 has absorption lines at this wavelength that 
are well suited to DIAL in absorption line strength, insensitivity to temperature, and freedom from interference by other 
gases. 

A holmium laser material, with an activator of thulium, is used to generate the 2-tim wavelength. The host material of 
yttrium lithium fluoride (YLF) has proven robust under strong pumping conditions with a low level of thermal lensing. 
Tests have shown as much as 600 mJ with a single frequency spectrum using Ho:Tm:YLF. 9 However, recent research 
has proven that the host material of lutetium lithium fluoride (LuLiF) can produce more energy that YLF for the same 
pumping configuration. 10 Validar is now being converted to operation with Ho:Tm:LuLiF, but the atmospheric 
measurements shown in Section 3 were taken using Ho:Tm:YLF. The new configuration with Ho:Tm:LuLiF should 
extend the range capability of the instrument with the availability of higher pulse energy. 



Reported Configuration* 

Upgraded Configuration 

Laser material: 

Ho:Tm:YLF 

Ho:Tm:LuLiF 

Pulse energy: 

75 mJ 

95 mJ 

Pulse width: 

150 ns 

140 ns 

Pulse repetition rate: 

5 - 10 Hz 

5 - 10 Hz 

Spectrum: 

single frequency 

single frequency 

Wavelength: 

2051.0 nm 

2053.5 nm 

Wavelength tuning: 

3-4 nm including C0 2 lines 

3-4 nm including C0 2 lines 

Beam quality: 

<1.3 times diffraction limit 

<1.3 times diffraction limit 

Long term (one hour) 



wavelength stability: 

< 16.5 MHz (for C0 2 DIAL application) 

not yet tested 

Wavelength accuracy: 

<15 MHz (for C0 2 DIAL application) 

not yet tested 

Detector 

InGaAs in heterodyne configuration 

InGaAs in heterodyne config. 

Telescope aperture: 

4 inches 

6 inches 

Scanner: 

8.5 inch aperture. 

8.5 inch aperture. 


full hemispherical coverage 

full hemispherical coverage 

Tablet: Specifications of Validar. 

*Data in Section 3 was taken with the Reported Configuration. 


The Ho:Tm:LuLiF rod is contained in a pump head using conductively cooled laser diode arrays. Water cooling is used 
for the heat sink of the diode arrays, but the laser diodes are not directly cooled with water. Six AlGaAs diode arrays 
side pump the laser rod; each array has six bars that can provide up to 600 mJ near 792 nm in a 1-ms pulse. Three sets 
of two diode arrays, side by side, are arranged 120° apart around the laser rod circumference. The laser rod has a 
diameter of 4 mm and is encased in a 6-mm outer-diameter fused-silica glass tube for water cooling of the laser rod; the 
coolant temperature is kept at 15 C. The 20-mm long laser rod is doped with 6% Tm and 0.4% Ho. The resonator is in 
a ring configuration with a total length of 2.0 m, with an acousto-optic Q-switch to spoil the laser cavity. 

A continuous wave (CW) Ho:Tm:YLF laser serves as both an injection seed source and a local oscillator. This laser 
was originally built by Coherent Technologies, Inc. for the Space Readiness Coherent Lidar Experiment (SPARCLE). 11 
Approximately 10 mW of this CW laser is split off and focused into a polarization maintaining optical fiber for the local 
oscillator channel. The rest of the CW energy is used for injection seeding after being frequency shifted 105 MHz by an 
acousto-optic modulator. Such an offset sets an intermediate frequency between the local oscillator and pulsed laser 
output. A series of three Faraday isolators are used between the CW laser and the pulsed laser in order to prevent 
damage to the CW laser from a counter-propagating laser pulse and to isolate the local oscillator channel from 
influences by the laser pulse and acousto-optic modulator. 






To ensure a single-frequency pulsed output by the injection seeding process, the pulse laser cavity is actively matched to 
the frequency of the injection seed by a ramp-and-fire technique. 12 A piezo-electric translator (PZT) mounted to a 
mirror of the pulsed laser is ramped by voltage, as shown in Figure 2. An InGaAs photodiode circuit views the 
resonance of the injection seed with the pulsed laser cavity, and resonance peaks occur when the pulsed laser cavity 
length matches the injection seed frequency. A circuit senses the resonance and fires the Q-switch when this cavity 
matching condition is determined. This results is a single-frequency output spectrum from the pulsed laser. Successful 
seeding can be determined by a smooth shape to the laser pulsed, as shown in Figure 2, and a doubling of output power 
as the ring resonator is forced to run in a single direction. 
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Figure 2: Signals involved in resonance matching. The lower 
trace is a ramp voltage applied to a PZT on a pulsed laser cavity 
mirror. Resonance occurs (upper trace) when the seed laser 
matches a mode of the pulsed laser. 


Figure 3: Single-frequency output pulse viewed in 
time domain (upper trace) and in frequency domain 
(lower trace). No mode beating is observed. 


Control of the wavelength of the CW laser has been shown to subsequently control the pulsed laser. For C0 2 
measurement applications, knowledge and control of the wavelength is critical. DIAL uses two wavelengths: one 
matched to an absorption line of C0 2 and one tuned off of the absorption peak. A technique has been developed for 
stabilizing CW 2 -pm lasers to absorption peaks by a wavelength modulation spectroscopic technique. 13 Modulation of 
the laser wavelength, and tuning control for feedback stabilization, is provided by a PZT mounted to the output coupler 
of the CW laser. Experiments showed that the laser’s frequency can be held to within a few percent of the absorption 
linewidth. Further research verified that the pulsed laser is closely locked to the seed laser wavelength so that the 
pulsed laser wavelength is stabilized to the absorption line as the seed laser is stabilized. 14 DIAL measurements also 
require that the wavelength on and off the absorption line be alternated. In the data of Section 3 and in a more complete 
report in preparation, we alternate wavelength at the pulse repetition rate of the laser by a square wave applied to the 
CW laser PZT. 15 Another possibility is to use a pulse format in which the laser is fired twice during a pumping cycle, 
giving a pulse separation of several hundred microseconds. 16 One pulse of this pair can be tuned on line and the other 
offline. To accomplish rapid wavelength switching in under 100 ps, we have used an electro-optic switch to select one 
beam or another generated by two different CW lasers. 


2.2. Receiver 


The pulsed laser output is transmitted to the atmosphere via a 6-inch diameter off-axis paraboloid telescope. After 
expansion by the telescope, the ouput laser beam can be pointed or scanned by a scanner mounted through the roof of 
the trailer. The outgoing pulse and atmospheric backscatter are separated by a polarization relationship imposed by the 
combination of a quarter-wave plate and polarizing beam splitter. 

Heterodyne detection is provided by InGaAs photodiodes in a dual-balanced configuration as diagrammed in Figure 4. 
The dual-balanced design is used here as it allows collection of all of the atmospheric backscatter and local oscillator 
through a 50/50 evanescent wave coupler. A more typical design would be to use a single photodiode with a coupler 
taking 90% of the atmospheric backscatter and 10% of the local oscillator. The dual-balance configuration has an 
additional benefit of canceling amplitude noise that may occur in the local oscillator. A photodetector circuit is being 
used that combines the photodiodes and amplifiers (in die form) integrated together in the same microelectronic 
package. Such packaging reduces parasitic capacitance, thus improving performance where a wide bandwidth is 
needed. A factor of 2.2 reduction in parasitic capacitance has been shown over a previous design using discrete 
components integrated on a printed circuit board. The circuit used in Validar has a 200 MHz bandwidth, but circuits 
have been developed with 2 GHz bandwidth. 
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Signal 


Figure 4: Diagram of integrated photodetector circuit. The components in the shaded region are integrated in a 
monolithic package. Polarization maintaining fibers deliver light to the InGaAs photodiodes. 


2.3. Data acquisition and signal processing 


Processing of the lidar data is accomplished with the system diagrammed in Figure 5. The heterodyne signal is 
conditioned by an analog front end that includes an anti-aliasing filter and two different level of amplification switched 
in to bring the outgoing scatter pulse and atmospheric return within the dynamic range of a digitizer. Sampling by the 
digitizer is at 500 Ms/s and 8 bit resolution. An example of a single digitized heterodyne signal is shown in Figure 6 in 
time domain. 

Digitized data flows to an array of SHARC digital signal processors (DSP) for real-time reduction of the heterodyne 
signal to wind and power data products. The DSP, digitizer, and control software are hosted by a single computer over 
a PCI backplane. Control software includes a program in which the user can set such parameters as number of sample 
points, range bin size, number of pulses to be averaged, scanner angles, and display preferences. An outline of the 
processing steps for an acquisition along a single line of sight is given in the following: 


Data Acquisition and Processing System 



Figure 5: Layout of the Validar electronics and real-time processing system. 


2.3.1. Step (A): acquire 

An individual pulse is digitized starting from a trigger generated by a monitor photodetector. A user-specified number 
of samples is taken and separated into range bins with an example setup being 50,000 samples broken into 512 sample 
bins. The first range bin, however, can be of a longer length of 1024 samples to isolate the outgoing pulse scattered 
from the lidar optics. The scatter pulse is also called a monitor pulse because it is used to monitor for frequency that 
may occur in the injection seeding process. Range bin information is converted to range in meters. A number of pulses 
will be acquired for averaging as specified by the user. 


2.3.2. Step (B): FFT 

The power spectrum of the data for each range bin is calculated by a fast Fourier transform (FFT). A Flamming window 
may be applied as an option. 


2.3.3. Step (C): monitor frequency determination 

The maximum magnitude frequency is estimated in the first range bin containing the monitor pulse. This estimated 
frequency is compared to a user-specified frequency range to check if the laser pulse was properly injection seeded. An 
unseeded pulse, which occurs for less than 1% of the laser pulses, will have a maximum frequency estimate at a 
frequency corresponding to the free spectral range of the pulsed laser cavity since multiple modes are beating. Such a 
frequency is well away from the intermediate frequency and is easily filtered by this step. If a bad pulse is detected it is 
discarded from the set so that it will not be included in averaging. 


2.3.4. Step (D): shift spectra 









Taking the first pulse as the reference, provided it passed the test of the above step, the following pulses are shifted in 
frequency to remove the effects of pulse-to-pulse frequency jitter. The amount of shifting applied to each pulse is 
determined from Step (C). 



Figure 6: A single pulse return. 


2.3.5. Step (E): averaging 

The jitter-corrected spectra for each range bin are averaged among all the pulses acquired. 

2.3.6. Steo (F): Doppler shit estimation 

The maximum power frequency is determined for each range bin and subtracted from the zero Doppler reference 
frequency. This result is used to determine the Doppler shift. 


2.3.7. Step (G): power estimation 

The signal power in each range bin is calculated. This information indicates the level of aerosol concentration and the 
presence of clouds. Power is also used in DIAL measurements in which the atmospheric backscatter is compared at two 
wavelengths to infer CCL concentration. Signal processing for coherent DIAL measurements is further described in a 
report currently in preparation. 15 


Each of the above steps may be repeated in different viewing angles. A typical measurement scenario is to steer the 
laser beam in three directions — two orthogonal azimuths for finding the horizontal wind field and zenith viewing for 
finding vertical wind motion. 


3. ATMOSPHERIC MEASUREMENTS 

Measurements were made of wind and C0 2 during the Spring and Summer of 2003 using a Ho:Tm:YLF laser with 75 
mJ output energy. Upgrades to the lidar, reflected in Table 1, are now being installed and should offer better 
performance in range capability. 


3.1. Wind measurements 

Vertical wind profiles are shown in Figure 7, with Figure 7a being the case of a cloudy day and Figure 7b being the case 
of a clear sky. In both these data sets the lidar is pointed toward zenith and wind data is continuously sampled over 
several minutes with 20 pulses grouped together for averaging and processing. The upper plot in each figure is the wind 
speed versus altitude and the lower plot is signal power versus altitude. Figure 7a shows a constant layer of cloud at 8 
km altitude and the passing of occasional cumulus clouds. The wind speed under these cumulus clouds tends to be 
upward at 1 m/s, showing the convection of air into the clouds. Aerosol structure can be observed in the signal power 
with the transition at around 2 km altitude from the atmospheric boundary layer (ABL) to the free troposphere marked 
by sharp decrease in signal by over an order of magnitude. As much as two orders of magnitude difference in 
backscatter between the boundary layer and free troposphere has been observed in other data sets. Velocity 
measurements in Figure 8a are made consistently up to 6 km before signal starts to fade into noise. Good measurements 
pick back up again from clouds starting at around 7.5 km. 

The example of Figure 7b was taken on a day when no clouds could be observed. Flow ever, Validar revealed the 
presence of sub-visible cirrus at 12.5 km altitude. Again, a sharp transition in backscatter is seen between the 
atmospheric boundary layer and free troposphere. Good velocity measurements in this case are made from aerosols up 
to 5 km altitude. The vertical wind speed is persistently near zero, as opposed to the day of convective activity of 
Figure 7a. 

A horizontal wind profde can be measured by viewing two orthogonal azimuths at the same elevation, determining 
speed along each azimuth, and vector summing the components as a function of altitude. In this approach the vertical 
wind component is assumed to be zero. An example is shown in Figure 8, taken with a 45 degree elevation angle. The 
wind vector is graphed in two parts as speed and direction. Aerosol backscatter gives good velocity measurements up to 
5.5 km. Spurious signals occur at 5.5 to 6 km from a loss of usable signal-to-noise ratio, such that the frequency 
estimator returns gibberish. A cloud at 6 km altitude, though, allows good velocity measurements to be made. 
Differences can be discerned from this data in the wind within the ABL and above the ABL, with the transition between 
the two occurring at 2 km altitude. In the ABL the wind shows a constant speed with changing direction, while the free 
troposphere shows an increasing speed with a constant direction. Increasing speed with altitude is to be expected as 
friction with the ground has less influence. 

The data acquisition and control software, described in Section 2, can automate the collection of data similar to Figures 
7 and 8. A typical measurement scenario is to view three different angles with 20 pulses averaged at each angle. At 5 
Hz pulse repetition rate such a measurement set requires 15-20 seconds to complete depending on the extent of data 
archiving selected. This measurement set, which provides horizontal and vertical wind profiles, can then be 
automatically repeated. 
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Figure 7: Vertical wind speed time history for (a) cloudy conditions and (b) clear sky. Velocity is shown 
in the upper graph of each data set. Power of the atmospheric return is the lower graph. 
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Figure 8: Horizontal wind profile. 


3.2. C0 2 measurements 

The configuration of Validar for C0 2 measurement is similar to that for wind measurements, with the addition of 
precise wavelength control and alternation among laser pulses between wavelengths on and off an absorption peak. 
Signal processing aspects are also different for DIAL in that backscatter power measurements are required rather than 
spectral shift information used for wind measurement. A detailed report of recent implementation of these technologies 
is in preparation, and only a brief overview is given here. 15 

One challenge in precise power measurement from a heterodyne signal is that the power return from a particular range 
depends on the phase of the atmospheric backscatter relative to the local oscillator. This effect, called speckle, creates a 
power fluctuation from pulse to pulse. Averaging of many pulses is required to reduce the effect of speckle — 1000 




pulses were averaged in our experiment at each of two wavelengths. Since the two wavelengths are alternated, a single 
DIAL measurement requires 2000 consecutive pulses. Processing of the DIAL data separates the 2000 pulses into pairs 
with, for example, the odd numbered pulses being the shots taken with wavelength on line and the even numbered 
pulses off line. Strict consecutive acquisition must be observed over these 2000 pulses, without skipping a pulse that 
woidd ruin the alternated wavelength sequence. 

Processing of the DIAL data uses the same steps for wind data described in Section 2, with the power estimation from 
each range bin being the information of interest. The best approach found for power measurement involves integrating 
frequency domain spectra over a threshold set just above the noise floor. Figure 9 is an example of this power 
measurement. The on-line case falls off at a faster rate with range than the off-line case as a result of strong 
atmospheric absorption on line. Comparison of the rate of attenuation of the on-line data and off-line data gives a 
measure of CCL concentration. Extensive tests were run to assess the precision in measuring CCL and to show that a 
diurnal variation of C0 2 could be detected. 15 



Figure 9: Power in the atmospheric return for on-line and off-line wavelengths. 1000 pulses were averaged for 
each wavelength. 
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